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ABSTRACT 

We consider a very large class of hierarchies of zero-curvature equations constructed from 
afhne Kac-Moody algebras Q. We argue that one of the basic ingredients for the appearance 
of soliton solutions in such theories is the existence of "vacuum solutions" corresponding 
to Lax operators lying in some abelian (up to central term) subalgebra of Q. Using the 
dressing transformation procedure we construct the solutions in the orbit of those vacuum 
solutions, and conjecture that the soliton solutions correspond to some special points in 
those orbits. The generalized tau-function for those hierarchies are defined for integrable 
highest weight representations of G, and it applies for any level of the representation and it 
is independent of its realization. We illustrate our methods with the recently proposed non 
abelian Toda models coupled to matter fields. A very special class of such theories possess 
a U(l) Noether charge that, under a suitable gauge fixing of the conformal symmetry, is 
proportional to a topological charge. That leads to a mechanism that confines the matter 
fields inside the solitons. 



1 Two talks presented by the authors at the "International Workshop on Selected 
Topics of Theoretical and Modern Mathematical Physics - SIMI/96", Tbilisi, Georgia, 
September/96. 



1 Introduction 



The study of soliton solutions of non linear differential equations has been 
developed considerably in the last decades using (apparently) quite diverse 
methods. In spite of the great variety of types of equations considered, some 
basic features seem to be common to a large class of them. Attempts to 
unify such various aspects of basic classical Soliton Theory is clearly very 
important and may lead to new insights into the role of solitons in Physics 
and Mathematics. 

Among the several methods of constructing solutions for non linear differ- 
ential equations we have the Hirota method M, the dressing transformation 
procedure [0, |[ |, [5|, ||, Backlund transformations j7|, ||, the inverse scat- 
tering method ||, the Leznov-Saveliev algebraic construction [ 10 1 and the 
tau- function approach |TT], |3|. Each one of these methods have their own 
advantages, and the choice of one or the other depends on the particular 
problem and model one wants to address. However, soliton solutions are 
special and the theories presenting them must possess some common struc- 
tures. By solitons here we mean a solution localized in space that travels 
without dispersion, and that keeps its form when scattered by other soliton, 
suffering just a shift in its position with respect to the one it would have if 
not for the scattering. In fact, in some theories like the Toda models which 
possess two dimensional Lorentz invariance, the solitons can have a particle 
interpretation and there are indications for the existence of a duality between 
the one-soliton solutions and the fundamental particles. In four dimensions, 
a similar duality, generalizing the electromagnetic duality, seems to exist 
between monopoles and gauge particles fill . 

Pratically all theories in one and two space time dimensions, presenting 
soliton solutions, have a representation in terms of a zero curvature condi- 
tion or Lax-Zakharov-Shabat equation. In addition, the corresponding Lax 
operators lie in some infinite dimensional Lie algebra. In fact, we can say 
that basically almost all known soliton equations are related to Kac-Moody 
algebras |H| [HJ. 

In this paper we try to unify some aspects concerning soliton solutions. 
We argue that a basic ingredient for the appearence of soliton solutions is that 
there must exist one or several solutions, which we call "vacuum solutions" , 
such that the Lax operators, when evaluated on them, should lie in some 
abelian subalgebra (up to central term) of the Kac-Moody algebra associated 
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to the model. Such subalgebra can be written as an algebra of oscillators bi 



[bi,bj] =i PiC 5 i+jfi (1) 

with fa being some complex numbers and C the central element of the Kac- 
Moody algebra Q. In several cases, but not in all, ([[]) constitutes a Heisenberg 



subalgebra |15j of Q. In addition, we argue that the components of the Lax 
operators in the direction of the 6j's should be constant for the vacuum 
solutions. 

Using the dressing transformation method one can then construct, out of 
a given vacuum solution, an orbit of solutions parametrized by a constant 
group element p of the Kac-Moody group. We conjecture that the soliton 
solutions of those theories lie on such orbits, and that they correspond to 
points of the orbits where p is the product of exponentials of eigenvectors of 
the constant elements defined by the Lax operators evaluated on the corre- 
sponding vacuum solution. Such observations provide not only a powerful 
and elegant method of constructing soliton solutions, but also allow us to 
connect and generalize several results known in the literature jl6| [17|]. The 
so called solitonic specialization of the Leznov-Saveliev solution proposed in 



lg, [L9[ |20[ in the context of Toda type models, can then be connected to the 
dressing method applied to such vacuum solutions. In addition, we believe 
that in several cases there is a connection with the Backlund transformation 
method. 

We can also connect that observation with tau-functions. We define tau- 
functions as states in integrable highest weight representations of the Kac- 
Moody algebra, lying in the orbit, under the action of the group elements 
performing the dressing transformations, of the highest weight state. Our 
definition is independent of the level of the representation and also on the way 
it is realised, constituting a generalization of the previous definitions of tau- 



functions for level one vetex operator representations ||TTJ. The connection 
with the Hirota method is then made by realizing the Hirota's tau-function 
as projections of the tau-function on some suitable states of representation. 
The truncation of the Hirota's expansion is then understood in terms of the 
nilpotency of some operators in those integrable representations. 

We illustrate our methods with the recently proposed non abelian Toda 



models coupled to matter fields [ 21fl . In fact, we discuss a very special class 



of such theories where the solitons play a crucial role. These models possess 
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a U(l) Noether charge that, under a suitable gauge fixing of the conformal 
symmetry, is proportional to a topological charge. That leads to a mechanism 
that confines the matter fields inside the solitons. 

The paper contains the material presented by the authors in their talks at 
the "International Workshop on Selected Topics of Theoretical and Modern 
Mathematical Physics - SIMI/96", held in Tbilisi, Georgia (September/96), 



and summarizes results of references [22, 21, 23|. It is organized as follows 



In section 2 we introduce the type of hierarchies of soliton equations we 
shall consider, and discuss their vacuum solutions. In section 3 we define 
the dressing transformations and in section 4 we construct the soliton solu- 
tions. Section 5 introduces the tau-functions and discuss how they connect 
to previous definitions and to the Hirota's tau functions. The non abelian 
Toda models are introduced in section 6, their soliton solutions constructed 
in section 7, and their properties discussed in sections 8 and 9. In section 
10 we present in great detail, the example of the model associated to the 
principal gradation of sZ(2)W. 

2 Hierarchies and vacuum solutions 

Non-linear integrable hierarchies of equations are most conveniently discussed 
by associating them with a system of first-order differential equations 

C N * = , (2) 

where are Lax operators of the form 

and the variables are the various "times" of the hierarchy. Then, the 
equivalent zero-curvature formulation is obtained through the integrability 
conditions of the associated linear problem (H), 

[jC n , C M ) = . (4) 

An equivalent way to express the relation between the solutions of the zero- 
curvature equations and of the associated linear problem is 
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The class of integrable hierarchies of zero-curvature equations that will be 
studied here is constructed from graded Kac-Moody algebras in the following 
way. Consider a complex affine Kac-Moody algebra Q = Q +€D , associated 
to a simple finite Lie algebra Q of rank r, and an integer gradation of its 
derived algebra Q labelled by a vector s = (s , si, . . . , s r ) of r+1 non-negative 
co-prime integers such that 

£=0&(s) and [&(s),^-(s)] cg i+J (s). (6) 



According to [24] , integral gradations of £ are labelled by a set of co-prime 



integers s = (s , si, . . . s r ), and the grading operators are given by 



Q s = H s + N s D--^-Ti(H s ) 2 C. (7) 



2iV s 



Here 



H s = Y,*aK-H°, N s = J2^t, ^ = £m^a a , mj = l; (8) 

a=l i=0 a=l 

H° is an element of the Cartan subalgebra of Q\ a a , a = 1,2, ...r, are 
its simple roots; ip is its maximal root; mf the integers in expansion ip = 
E„ =1 mja a ; and are the fundamental co-weights satisfying the relation 

We have in mind basically two types of integrable systems. The first 
one corresponds to the Generalized Drinfel'd-Sokolov Hierarchies considered 
in [|14|| , and P5| , which are generalizations of the KdV type hierarchies studied 



in In particular, and using the parlance of the original references, we 
will be interested in the generalized mKdV hierarchies, whose construction 
can be summarised as follows (see and, especially, ^ZE \ for details). Given 



an integer gradation s of Q and a semisimple constant element E\ of grade I 
with respect to s, one defines the Lax operator 

L = d x + Ei + A , (9) 



where the components of A are the fields of the hierarchy. In |25[, it was 
shown that the component of A along the central term of Q should not 
be considered as an actual degree of freedom of the hierarchy. This is the 
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reason why these hierarchies can be equivalently formulated both in terms 
of affine Kac-Moody algebras or of the corresponding loop algebras. They 
are functions of x and of the other times of the hierarchy taking values in 
the subspaces of Q with grades ranging from to I — 1. For each element 
in the centre of Ker(ad E{) with positive s-grade N, one constructs a local 
functional of those fields, B N , whose components take values in the subspaces 
^o(s), . . . ,Q N (s). Then, B N defines the flow equation 

M- = [Bn,L], (10) 
ot N 

and the resulting Lax operators = d/dt^ — B^ commute among them- 



selves [14 1 



The second type of integrable systems corresponds to the non-abelian 
affine Toda theories |^6[ |23|, |21], pTfl , and a very general class of these models 



will be described in section ^. 

An important common feature of all those hierarchies is that they possess 
trivial solutions which will be called "vacuum solutions". These particular 
solutions are singled out by the condition that the Lax operators evaluated 
on them lie on some abelian subalgebra of Q , up to central terms. Then, the 
dressing transformation method can be used to generate an orbit of solutions 
out of each "vacuum". Moreover, it is generally conjectured that multi- 
soliton solutions lie in the resulting orbits. As a bonus, the fact that we 
only consider the particular subset of solutions connected with a generic 
vacuum allows one to perform the calculations in a very general way and, 
consequently, our results apply to a much broader class of hierarchies. 

For a given choice of the Kac-Moody algebra Q and the gradation s, 
let us consider Lax operators of the form (^) where the potentials can be 
decomposed as 

N+ 

A N = ^2 A N,i , where A Nji G Gi(s) (11) 

i—N— 

N- and iV + are non-positive and non-negative integers, respectively, and 
the times tpj are labelled by (positive or negative) integer numbers. The 
particular form of these potentials will be constrained only by the condition 
that the corresponding hierarchy admits vacuum solutions where they take 



5 



the form 

N- 

I ( vac ) 



^ aC) = E + Mt) c = e N + f N (t) C . (12) 



In this equation, C is the central element of Q, and 6j G Gi(s) are the gener- 
ators of a subalgebra s of Q defined by 

s = {bi e Qi(s) , i e E c1L \ [b u b,] = i^C 5 i+jfi } , (13) 

where are arbitrary (vanishing or non-vanishing) complex numbers such 
that /3_j = and i£ is some set of integers numbers. Moreover, also 
arbitrary numbers, and /jv(t) are (D-functions of the times i/v that satisfy the 
equations 

df N (t) df M (t) 



dt M dt N 



These vacuum potentials correspond to the solution of the associated 
linear problem given by the group element (|5|) 

^ {vac) = exp e A + T(*) CJ (15) 
where the numeric function ^(t) is a solution of the equations 

°' Af] -f N (t) + lj2ip i( f N c^t M . (16) 
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3 Dressing Transformations 

In terms of the associated linear problem, one can define an important set 
of transformations called "dressing transformations", which take known so- 
lutions of the hierarchy to new solutions. Regarding the structure of the 
integrable hierarchies, these transformations have a deep meaning and, in 
fact, the group of dressing transformations can be viewed as the classical 
precursor of the quantum group symmetries |J. Denote by G _(s), G +(s), 
and G (s) the subgroups of the Kac- Moody group G formed by exponen- 
tiating the subalgebras £< (s) = © i<0 £i(s), £>o(s) = ® i>0 Gi{s), and G (s), 
respectively. According to Wilson [|l(| [HJ, the dressing transformations can 
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be described in the following way. Consider a solution \1/ of the linear prob- 
lem (H), and let p = p_ po P+ be a constant element in the "big cell" of G , 
i.e., in the subset G _(s) G Q (s) G + (s) of G , such that 

* p = p ^-i) <0 (^ p $-i) ($ p *-i) >0 . (17) 

Notice that these conditions are equivalent to say that both p and ^ p ^l/ -1 
admit a generalized Gauss decomposition with respect to the gradation s . 
Define 

w = e {0) {(^ p^- 1 )^]- 1 ^p 

= eV 0) (v[/pvI/- 1 ) >0 vD (i 8 ) 

where 

e 10 ' ef'^fpr 1 )^ (19) 

Then, \l/ p is another solution of the linear problem. In order to prove it, 
introduce the notation g± = p ^l/ -1 )^ and 9^ = d/dtjv, and consider 

+ eL ^:: 1 ^**" 1 )'/-©- 5 " 1 

+ ei 0) ff+ (^* *- 1 )^ 1 ei 0> " 1 (20) 



Then, the first identity implies that 3n^ p ^ p 1 G ©j<7v + ^i( s ); an d the second 
that 9iv* p G 0i>7v_ ft(s). Consequently 

^ = 757- ^ e © ^( s ) > ( 21 ) 
ot N i=N _ 

and, taking into account (|TTJ), it is a solution of the hierarchy of zero- 
curvature equations. If the fields of the hierarchy are such that A^^ does not 
span the whole subspace Qi(s) then we have to impose further constraints on 
the group elements performing the dressing transformation. 
For any p lying in the big cell of G , the transformation 

V p : h+ W , or A N ^A P N , (22) 

is called a dressing transformation. 
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4 Solitons out of vacuum solutions 



We now consider the orbit of the vacuum solution (|15|) under the group of 
dressing transformations. For any element p of the big cell of G , let us define 

6 > = ^(vac) p ^(vac)-^ @< = ^(vac^vac)-^ ~ X ( 2 3) 

and 

e + = efe>, e_ = 0(° ) 0<; (24) 

where 0+^ and 0^ are the same as in fll9|), but with ^ replaced by \l/( vac ) 
defined in flloT). 



Now, the orbit of the vacuum solution (FJ) can be easily constructed 
using Eqs. ([Tj|) and (|2l"l). Under the dressing transformation generated by p, 



or, equivalently, A^ ac ' ) becomes 

^-/iv(t)c = 0_ £at 0_ _1 + <9jV©- 0_ 1 G &(s) 

i<7V+ 

= + e 7V e + - 1 + 0^8+ 6+" 1 G &(s), (26) 

i>N- 



Eqs. (p4|) , ( p3|) and ( p6|) summarize the outcome of the dressing transfor- 
mation method, which, starting with some vacuum solution (fl2|) , associates 
a solution of the zero-curvature equations (|J) to each constant element p in 
the big cell of G . The construction of this solution involves two steps. First, 
the Eqs. ( PB"D can be understood as a local change of variables between the 
components of the potential An and some components of the group elements 
0+ and 0_. 

The second step consists in obtaining the value of the required com- 
ponents of + and 0_ from Eqs. (p4|) and (f23f). This is usually done by 
considering matrix elements of the form 

(p | 0+ | fi') = (p | e ^" ENtN pe~^N ENtN | //) , (27) 

where | p) and | p') are vectors in a given representation of Q. The ap- 
propriate set of vectors is specified by the condition that all the required 
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components of G_ and + can be expressed in terms of the resulting matrix 
elements. It will be shown below that the required matrix elements, consid- 
ered as functions of the group element p, constitute the generalization of the 
Hirota's tau-functions for these hierarchies. Moreover, Eq. ( p7[) is the ana- 
logue of the, so called, solitonic specialization of the Leznov-Saveliev solution 
proposed in [19], 23], for the affine (abelian and non-abelian) Toda 



theories. 

Consider now the common eigenvectors of the adjoint action of the ejy's 
that specify the vacuum solution ([12|) . Then, the important class of multi- 
soliton solutions is conjectured to correspond to group elements p which are 
the product of exponentials of eigenvectors 

p = e Fl e F * ... e F " , [e N , F k \ = uff F k , k = 1, 2, . . . n . (28) 

In this case, the dependence of the solution upon the times tw can be made 
quite explicit 



(k) 

(p | e: 1 e + | //) = (p \ J[ exp(e£>* tn Fk) \ ^ . (29 ) 

k=l 



We emphasize that not all solutions of the type (|29|) are soliton solutions, 
but we conjecture that the soliton and multi-soliton solutions are among 
them. The conjecture that multi-soliton solutions are associated with group 



elements of the form (j28|) naturally follows from the well known properties 
of the multi-soliton solutions of affine Toda equations and of hierarchies of 
the KdV type, and, in the sine-Gordon theory, it has been explicitly checked 
in Ref. ||. Actually, in all these cases, the multi-soliton solutions are ob- 
tained in terms of representations of the "vertex operator" type where the 
corresponding eigenvectors are nilpotent. Then, for each eigenvector F k there 
exists a positive integer number such that (i 7 ^)" 1 7^ only if m < m^. 
This remarkable property simplifies the form of ( ppj ) because it implies that 
e Fk = 1 + .Ffc + ■ • ■ + (Fk) mk / m^l , which provides a group-theoretical 
justification of Hirota's method. 

An interesting feature of the dressing transformations method is the pos- 
sibility of relating the solutions of different integrable equations. Consider 
two different integrable hierarchies whose vacuum solutions are compatible, 
in the sense that the corresponding vacuum Lax operators commute. Then, 
one can consider the original integrable equations as the restriction of a 
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larger hierarchy of equations. Consequently, the solutions obtained through 
the group of dressing transformations can also be understood in terms of the 
solutions of the larger hierarchy, which implies certain relations among them, 
(see section 4 of fl23| for more details). 



5 The tau-functions 

According to the discussion in the previous section, the orbits generated 
by the group of dressing transformations acting on some vacuum provide 
solutions of certain integrable hierarchies of equations. Making contact with 
the method of Hirota, the generalized "tau-functions" that will be defined in 
this section constitute a new set of variables to describe those solutions. One 
of the characteristic properties of these variables is that they substantially 
simplify the task of constructing multi-soliton solutions PB| . The group- 



theoretical interpretation of this property has already been pointed out in 
the previous section. Tau-functions are given by certain matrix elements in a 
appropriate representation of the Kac-Moody Group G . Moreover, the tau- 
functions corresponding to the multi-soliton solutions are expected to involve 
nilpotent elements of G , which is the origin of their remarkable simple form. 

The tau-function formulation of the Generalized Drinfel'd-Sokolov Hier- 
archies of [|14[] has already been worked out in [25], which, in fact, has largely 



inspired our approach. However, there are two important differences between 
our results and those of |25|] . Firstly, our approach applies to the affine Toda 
equations too, and, secondly, it does not rely upon the use of (level-one) 
vertex operator representations. 

At this point, it is worth recalling that the solutions constructed in sec- 
tions |3| and f| are completely representation-independent. In contrast, our 
definition of tau-functions makes use of a special class of representations of 
the Kac-Moody algebra Q called "integrable highest-weight" representations. 
The reason why these representations are called "integrable" is the following. 
For an infinite-dimensional representation, it is generally not possible to go 
from a representation of the algebra Q to a representation of the correspond- 
ing group G via the exponential map x ^ e x . However, the construction 
does work if, for instance, the formal power series terminates at a certain 
power of x, or if the representation space admits a basis of eigenvalues of x. 
These conditions, applied to the Chevalley generators of Q, single out this 
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special type of representations. 

The generalized tau-functions will be sets of matrix elements of the form 
indicated on the right-hand-side of ( P7| ) , considered as functions of the group 
element p. They are characterized by the condition that they allow one 
to parameterize all the components of G + and G_ required to specify the 
solutions (|26|) of the zero- curvature equations (|). As we have discussed 
before, the tau-functions corresponding to the multi-soliton solutions are 
expected to have a very simple form. However, in contrast with the original 
method of Hirota, we cannot ensure in general that the equations of the 
hierarchy become simpler in terms of this new set of variables. 

First, let us discuss the generalized Hirota tau-functions associated with 
the components of B. In equation (PT|), these components can be isolated by 
considering the vectors | /io) of an integrable highest-weight representation 
L(s) of Q which are annihilated by all the elements in (? >0 (s), i.e., T | /i ) = 
and (/i | T' = for all T e £>o( s ) an d T" G £/ <0 (s), respectively. Then, the 
corresponding tau-functions are defined as 

W*) = (K> I ^ (vac) p * (vac) ~ 1 1 /i ) 

= (// I e ^N £NtN pe~^* £NtN | /i ) , (30) 
and, in terms of them, equation Q27| ) becomes just 

(// | B~ l | /i ) = T m> ^(t) 

where we have denoted 

B- 1 = (\[/( vac ) p -qjt™^ 1 ^ 

By construction, ^o( s ) always contains the central element C of the Kac- 
Moody algebra, but it is always possible to split the contribution of the 
corresponding field in ( [31]) . Let s g ^ and consider the subalgebra of 
Q generated by the ef with i = 0, . . . , r but i ^ q, which is a semisimple 
finite Lie algebra of rank r (Q^ is always simple if q = 0). Then, £/o( s ) = 
(^o( s ) H Q^) ©CC and, correspondingly, B can be split as B = b exp(z/ c). 

2 Since the resulting relations between tau-functions and components of the j4at's will 
be considered as generic changes of variables, we will not generally indicate the intrinsic 
dependence of the tau-functions on the group element p. 
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(31) 



(32) 



Here, v is the field along C, and b is a function taking values in the semisimple 
finite Lie group G (q) whose Lie algebra is £? (s) n G {q) . Since K = £[ =0 fcV^ 
is the level of the representation L(s), Eq. ( ]5Tj) is equivalent to 



(/i'o | fl" 1 | /i ) = e"** (/i' | 6" 1 | no) = W£) (t) . (33) 

Moreover, it is always possible to introduce a tau-function for the field v. Let 
us consider the highest- weight vector | v q ) of the fundamental representation 
L(q), which is obviously annihilated by all the elements in Q^ q \ Therefore, 

(v q I I v 9 ) = e"^ v = r VqiVq (t) = rf (t) , (34) 



which leads to 



<// I 6" 1 | aio) = T^^v and „ = - In ^ . (35) 



q 



Finally, recall that the vectors | /xq) form a representation of the semisim- 
ple Lie group Go ■ Therefore, if L(s) is chosen such that this representation 
is faithful, Eq. (|35|) allows one to obtain all the components of b in terms 
of the generalized tau-functions r^y and . Notice that, in this case, 
the definition of generalized tau-functions coincide exactly with the quanti- 
ties involved in the solitonic specialization of the Leznov-Saveliev solution 



proposed in [20 



Let us now discuss the generalized tau-functions associated with the com- 
ponents of 0f . Consider the gradation s of Q involved in the definition of the 
integrable hierarchy. For each Si ^ 0, let us consider the highest-weight vec- 
tor of the fundamental representation L(i) and define the (right) tau-function 
vector 

| r p(t)) = # (vac) h ^( vac ) _1 | Vi ) 

Notice that | if(t)) is a vector in the representation L(i). Therefore, it has 
infinite components, and it will be shown soon that the role of the Hirota 
tau-functions will be played by a finite subset of them. Taking into account 
that | Vi) is annihilated by all the elements in g > o(s), equation ( P7| ) implies 

0< _1 B- 1 \ Vi ) =| T?(t)) , z = 0,...,r and s t ^ . (37) 
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The definition ( |36| ) is inspired by the tau-function approach of [IT , p5 



|23fl . However, in j25|1 , and |p3| ) the authors consider a unique tau-function 
| T s (t)) G L(s). In fact, one could equally consider different tau-functions 
| r s /(t)) associated with any integrable representation L(s') such that 
if, and only if, Sj 7^ 0. Since the highest weight state of L(s) is obtained by 
the tensor product of the | i>i)'s, as 

l«-> = ®{| «*>•"}■ (38) 

one sees that all these choices lead to the same results, but the one presented 
here is the most economical. 

Since, for any integrable representation, the grading operator (0) can be 
diagonalized acting on L(s), these tau-functions vectors can be decomposed 
as 

I if®) = £ I T« H \t)) , Qi I rf ^(t)) = -j I rf ^(t)) , (39) 
-ieZ<o 

where we have used that 0f G G <o(s) and 5 6 Go(s), and Qi indicates the 
derivation corresponding to the grading operator (|7|) with Sj = 5^%. More- 
over, the highest- weight vector is an eigenvector of the subalgebra &o(s) an d, 
consequently, of B. Therefore, 

K R(0) (t)> = B- 1 \ Vi ) = rP(t) \ Vi ), (40) 

where, t\ it) is a (D-function, not a vector of L(i), whose definition is f\ 

T?\t) = ( Vi I e X>"twp e-Etf*"** I Vi ) = T Vi>Vi (t) (41) 

(compare with Eq. (|33|)). Therefore, Eq. (^) becomes 

©r 1 I Vi ) = I 7?{t)) , (42) 



To compare with (|3l|), notice that | /io) =| fj) forms a one-dimensional representation 
of Go(s) and, consequently, ^^(i) vanishes unless | (i' ) =| «,-). Therefore, for non- 
abelian Go, the required tau-functions t^^' (t) have to involve the fundamental integrable 
representations L(j) corresponding to Sj = 0, in contrast with | r^it)) (see Eq. (|37 
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which is the generalization of the Eq. (5.1) of |25| for general integrable 



highest-weight representations of Q. Eq. (|42D allows one to express all the 
components of 05 in terms of the components of | T^{t)) for all % = 0, . . . , r 
with Si (for instance, by using the positive definite Hermitian form of 
L(i)). However, it is obvious that only a finite subset of them enter in the 
definition of the potentials A N through Eq. (p6|). 

In exactly the same way, one can introduce another set of "left" tau- 
function vectors through 



which leads to 



(r t L (t) | = ( Vi | ^ (vac) h ^ (vac) , (43) 

1 



(vi I ®l = tf(t) | , (44) 

and allows one to express all the components of 0> in terms of the compo- 
nents of (tj (t) | for all i — 0, . . . , r with 7^ 0. 

Summarising, the generalized Hirota tau-functions of these hierarchies 
consist of the subset of functions t M0! ^ and of components of | r/ 2 ) and 
(rf I required to parameterize all the components of the potentials in 
Eq. (|26|) . Then, for the multi-soliton solutions corresponding to the group 
element p specified in (|28|), their truncated power series expansion follows 
from the possible nilpotency of the eigenvectors F& in these representations. 
For instance, if n = 1 in and F[™ \ Po) — F™ I v i) = unless m < mi, 



then 



mi 1 

2 tt e * 04 I F* |/io> 



mi 



and 



mi 1 

\ Tl R (t)) = Tv ek ^ N UJNtN F i I v i) ■ ( 45 ) 

6 The example of the Non Abelian Toda 
Models 

Consider an untwisted affine Kac- Moody algebra Q endowed with an integral 
gradation Q = ® ne ~%_Gn (see ©,(0)). By an affine Kac- Moody algebra we 
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mean a loop algebra corresponding to a finite dimensional simple Lie algebra 
Q of rank r, extended by the center C and the derivation D. 

Let Ai be a two dimensional manifold with local coordinates x + and 
x_; Q be an affine Kac- Moody algebra corresponding to a finite dimensional 
complex simple Lie algebra Q with the Lie group G; A be a flat connection 
in the trivial holomorphic principal fibre bundle Ai x G i — ► Specify 
the connection in such a way that its (1, 0)-component takes values in the 
subspaces ® l n=0 G+ n , and (0, l)-component takes values in @ l n=0 G-n-> with 
I being a fixed positive integer. In other words, up to a relevant gauge 
tranformation, these components, satisfying the zero curvature condition 

d + A_ - 3-A+ + [A+ , AJ = 0, (46) 

are of the form 

A + = -B F + B- 1 , A_ = —d-B B l + F~ . (47) 

Here B is a mapping from .M to the Lie group Go with the Lie algebra Qq; 
F ± are mappings to 0^ =1 Q± n of the form 

F+ = £ z + F- = E_ /+ ^F-, (48) 

m=l m=l 

with E±i being some fixed elements of Q±i\ and F^, 1 < m < I — 1, take 
values in Q± m . 

Substituting the gauge potentials fl4"T|) into (pEBP , one gets the equations 
of motion 

i-i 



d+(d-BB- 1 ) = [E^ , B E t B- 1 ] + , B F+ B- 1 ] , (49) 

71=1 

Z— m— 1 

d_F+ = [Ek, B- 1 F l Z m B}+ l^+miB-'F-B], (50) 

71=1 

I — 771— 1 

«9 + F- = -[^.SF^Bl- 2 t F n + -' BF+B^).(51) 

71=1 

Since Q s , defined in (^), and C are in C/ , we parametrise £> as 

B = be r > Qa e v °, (52) 
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where b is a mapping to Go, the subgroup of Go generated by all elements of 
Go other than Q s and C. The fields r] and v correspond to the extension of 
the loop algebra, and, as we will show below, are responsible for making the 
system conformally invariant |28|, |29|]. Clearly, the order of the three factors 



in ( p2|) is irrelevant, since they commute. In addition, we will use a special 
basis for the generators of Go such that they are all orthogonal to Q s and 
C. From (|7|) one observes that the generators of Go are, besides C and Q s , 
the elements H®, a = 1, 2, ... r, of the Cartan subalgebra, and step operators 
E± a and E^l, such that YZ,=i s a^a ■ a = 0, and YZ,=i s a^a ' ft = N s - There 
can be no step operators E", with | n |> 1, as explained in appendix C of 



ref. | 2~3fl . Therefore, shifting the Cartan elements as 

H° a =H°-±- Tr (H.Hi) C = H° a - ±%-C, (53) 

one gets 

Tr (C 2 ) = Tr (c tf a °) = Tr (Q 2 S ) = Tr (Q s H°) = 0, Tr (Q s C) — N B , 
Tr (H° a SJ) = Tr (flj = 4a a ■ a b /a 2 a a 2 b = Vab , (54) 

for all a, b = 1 . . . , r. 

Here we have used H® = 2a a -H / a 2 , Tr {x ■ H° y ■ H°) = x-y, and Tr {C D) = 
1. For more detail of such a special basis, see appendix C of ref. p3[| . 
Substituting ( |52T) into the equations of motion (^9|)-(|5lD , one has 

i-\ 

d+fd^bb' 1 ) + d + d^uC = e lv [E^ l ,bE l b~ 1 ] + J2e nr '[F-,bF+b- 1 ](55) 

n=l 

l—m—X 

<9_F+ = ef-")"^,, 6-^^61+ £ e^[F+ +n , 6-^-6], (56) 



n=l 
Z— m— 1 



n=l 



d+d_r)Q s = 0, (58) 

where the last equation is a consenquence of the fact that D, and hence Q s , 
can not be obtained as the Lie bracket of any two elements of G ■ 

The structure of the vacuum of the system (|55|)-(|58|) is rather compli- 
cated. We will discuss some aspects of it below. However, there is a simple 
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condition that guarantees the existence of static (vacuum) solutions. If the 
elements E±i satisfy the relation 

[Ei , E_i] = (3C , where (3 = £ Tr {E l £_,), (59) 

then 

b = l, F± = 0, = , v = -f3x + x-, (60) 

is a (vacuum) solution of (|55|)-(|58|). 

Another possibility for vacuum solutions arises when E±i, I > 1, belong 



to a Heisenberg subalgebra of Q, see p3l 115 



[E M , E N ] = Tr (E M E_ M ) M S m+ n,o C, (61) 

where M,N belong to some (infinite) subset 2^ of the integer numbers 2. 
In such cases one has that 

6 = 1, 77 = 0, F± = c± E± M , F± = 0, m^TL E , v = -Qx + x_, (62) 

is a solution of (|55l) -(|58"l) with being constants, and 

= + ]T Tr (E M E_ M ) Mc + M c M . (63) 

M=\ 

Obviously, the system Q55])-(|5B|) may have many more vacuum solutions 
besides Q60|) and ( |6^) . However, the condition (|59|) guarantees the existence 
of at least one vacuum solution. Such a fact, as we will see below, favors the 
existence of soliton solutions. 

The models introduced above are completely characterised by the data 
{G, Qs, I, E±i}; and we have a quite large class of systems with physical prop- 
erties crucially depending on a choice of those data. 

Equations (^) - (|5*%D are invariant under the conformal transformation 

z+^f(x+), #(£_), (64) 

with / and g being analytic functions; and with the fields transforming as 

b{x + , xJ) — ► b(x + , xJ) = b(x + , x_) , (65) 

e -r?(a; + ,x_) _^ g-^(5+ _ ^ f'^ 1 (g') 1 ^ 1 e~ V ^ X+ ,X ~^ , (67) 

F+(x + ,x_) - F+(£ + ,£_) = (/')" 1+m// ^(a: + ,x_), (68) 
F~(x + , x_) -> F~(x+ , x-) = {g')' 1+m/l F~(x + , x_) , (69) 
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where the conformal weight 5, associated to e~ u , is arbitrary. 

Notice that the Lorentz transformation x± — > X Tl x± is obtained from 
Q5%D by taking f(x + ) = x + /\ and g(xJ) = Ax_. 

Equations (p5|)-([58|) are also invariant under the transformations 

b(x + ,x^) -> h L (x-) b (x + , x-) h R (x+) , (70) 
F+(x+,x_) -> ^(zh-) ^(s+j z-) h R (x+), (71) 
F m (x+,x_) -> /i L (x_) F" x_) /i L x , (72) 

where h L (x^) and /ir ( elements of subgroups 7Yq and of G 0) 

respectively, satisfying the conditions 

h R (x + ) E l h R 1 (x + ) = E l , hr L x {x„) E_ l h L {xJ) = E_ l . (73) 

The left and right gauge transformations commute, and so the gauge group 
is Hq ®TLq. Whenever 71$ and Hq have a set of common generators, we get 
an important subgroup of the gauge group, namely TCd = ^ ^-o- These 
are global gauge transformations, where the fields are transformed under 
conjugation (Ul = = hp = const.), 

b^h D bh D \ F±^h D F±h~ D \ (74) 

and E±i = hf ) E±ih] D 1 . We discuss the relevance of these transformations 
below. 



7 Soliton Solutions 

We now perform the dressing transformation, dicussed in sections |3| and f|, 
by taking as an initial configuration a vacuum solution of fl55| )-(|58|). As we 
have said, the model under consideration may have several type of vacuum 
solutions. However, here we will deal with the solutions of type ( |60"D or fl62|) . 
For the vacuum solutions the gauge potentials (pT7p become 

= -8+ , A {0) =£_+ Qx + C, (75) 

with 8± given by 

i-i 

£± = E ±l +J2 c%E ±N , and so [8+ , £-] = Q C; (76) 

N=l 



18 



where cfj and Vt were introduced in (p2|) and (|63|) , respectively. 
They can be written as 

Af = -d±^ iv£ic) ^( vac )"\ with ^ (vac) = e x+£+ e~ x - £ -. (77) 

The gauge potentials for the vacuum solution ( joTf ) are obtained from 
by taking = 0. In fact, they are connected by the gauge transformation 

A (0) = vj)(vac) A {0) |^ ^(vac)^ 1 _ ^(vac) ^(vac)^ 1 ? ^g) 

with 

^( vac ) = exp[o: + (S+ - E t )] exp[-a:_ (£_ - £_,)]. (79) 



However, in general, the vacuum solutions ((6^) and ( 62[) may not be con- 
nected by any dressing transformation, and, in such the existence of 
two elements of form (]23|), is not always possible. Consequently, one can have 
soliton solutions lying on different orbits under the dressing transformations. 

In order to perform the dressing procedure, we take Q75D as initial gauge 
potentials. Then, we obtain, under the dressing procedure, the solutions on 
the orbit of vacuum (|62[) , and for = those on the orbit of the vacuum 
fl60l) . From the structure of the dressing transformations and from the fact 
that the grading operator (|?]) is never the result of any commutation, since 
it contains D, it follows that the dressing transformations do not excite the 
field 77. Therefore, from ©, (H), © and we get 

bE l r 1 + b- 1 = 6± + Q ± d +®±) ®i\ (8°) 

m=l \ n=l I 



l-l 

- cL bb- 1 -(d„v + Qx+) C + + ]T F m = (81) 

771=1 

i-1 



e± [e-x + £ c - s_ m - e^cte^ q ± \ 
\ 777=1 / 



Note that in the above relations, the fields b, v and stand for the solutions 
on the orbit of the vacuum solution (|62|) . The procedure to construct the 
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solution requires to split the above equations into the eigensubspaces of the 
grading operator flj]). It is convenient to write 

6> = exp (]Tt (s) ) , 6< = exp [J2 t{ ~ s) J > where ti±s) e G±s- (82) 

\s>0 / \s>0 / 

The mappings t^ ±s \ for each choice of p, are determined from (1231) with 



\l/ being \]>( vac ) given in fl77|). Then, the components of ( j30|) and (82) in 
each eigensubspace, give an equation connecting the fields with t^ ±s \ Thus 
the solutions for the fields b, v and are determined from t^ ±s \ Such a 
procedure is rather cumbersome, but fortunately, one needs to know very 
few t( ±s )'s to get the solution. For instance, taking relations ( p0| ) and (|32]) 
for + (0_) with grade components and —I (I and 0), one gets 

ef = hi\x-), eL 0) = be {i/+Qx+x - )c h R (x + ), (83) 

with hi{x-) and /ir(x + ) defined in ([73|). 

From (|I7D, @, fl77|) and (||) it follows that 

0< _1 (/i L (x_)6e (l/+n:,;+x - )c /i i? ,(x + )) _1 e> = e x+£ +e- a: - £ -pe a: - £ -e- a: + £ +. 

(84) 

The space-time dependence of the r.h.s. of the above relation is given 
explicitly. One can extract the solutions out of (|8~||) by taking the expectation 
value of its both sides between suitable states of a given representation of Q, 
in a similar way to that one explained in section |j. 

The solitons solutions are obtained from ( |3~4"D by choosing the fixed group 
element p, characterising the dressing transformation, as the exponential of 
an eigenvector of £±, i.e. 

p = e v . (85) 

That is the solitonic specialization discussed in section |]. Indeed, if V satisfies 
the relations 

[£± ,V]=u± V, (86) 

then (Hp reads as 

exp (e x + u +~ x - w " v) = exp (e 1 ^'^ V) , (87) 
with 7 = uj + + UJ-, and v = (cj_ — uj + ) / (uj + + cj_), since x± — t ± x. 
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Therefore, for each eigenvector V, expression (R^) corresponds to a solu- 
tion that travels with a constant velocity v without dispersion. Depending 
upon the properties of V, as we will see below in the examples, such solutions 
correspond to one-soliton solutions. 

The multi-soliton solutions are obtained by taking p to be the product 
of several one-soliton p's, i.e., 

p = e Vl e Va e V3 ...e Vir , (88) 

with each Vi satisfying [£± , V*] = uj l ± Vi. 

Notice that, under the global gauge transformations ([71P, the gauge po- 



tentials (|47 ) are transformed as A± — > hp A± . Therefore, since the poten- 
cials are pure gauge, A^ = — c^TT -1 , one has T — > h^T, and consequently 
( |23|) implies 0+ — > hnQ+h^ 1 and B5 — > ho®- hp 1 . Hence, with solution 
( I81D corresponding to a fixed element p, a solution, obtained from that by a 
global gauge transformation (|74l) , is given by ( |34D with the replacement 

p^h D ph- D \ (89) 

if the condition hr,£± h^ 1 = £± is satisfied. For the solutions on the orbit 
of the vacuum (0), that is indeed true, since £± = E±i\ see (|76|) . For the 
solitonic case, one then obtains for each eigenvector V of £±, an orbit of 
equivalent one-soliton (or multi-soliton) solutions generated by hnVh^ 1 . 



8 Masses of fundamental particles and 
solitons 

As we have seen above, the system under consideration is conformally in- 
variant. Therefore, since we do not have a continuum mass spectrum, its 
fundamental particles have to be massless. However, such a symmetry can 
be spontaneoulsy broken by choosing a particular constant solution for the 
field 77, say rj = f] . The resulting theory is then massive. Representing the 
mapping B as B = expT, and considering only the linear field approxima- 
tion, i.e., the free part of the equations of motion fl49D-fl5lD, one gets 

d + d^T = [E^T]], (90) 

8 + d^F+ = —v v [E-i , [Ei , , (91) 
8 + d-F~ = —v v [E-i , [Ei , F~]] , (92) 
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where v v = e lTJ0 . 

Therefore, the masses of fundamental particles in such a theory are given 
by the eigenvalues of the operator [E_i , [Ei , *]] in the subspaces Q n , n = 
0,±1,±2,...±(/-1), i.e., 

[E-i , [E l ,X]] = XX. (93) 

Since <9 + <9_ = \{d 2 — d 2 ), we obtain the masses from the Klein-Gordon type 
equations (pOf) - ( |92"1) as 

m 2 x = A\v v . (94) 

That result constitute a generalization of the arguments used in the abelian 
and non abelian affine Toda models |3(| [23| . Of course, we are interested in 
those cases where the eigenvalues of the operator [E_i , [Ei , *]] are real and 
positive on the subspaces under consideration. That will be, in fact, one of 
the conditions we use to select the data {Q, Q s , I, E±i} for defining physical 
models through (fT^). 

Notice that the field e lr] plays the role of a Higgs field, since it not only 
spontaneously breaks the conformal symmetry, but also because its vacuum 
expectation value sets the mass scale of the theory. We have here the same 



mechanism as in non abelian affine Toda theories [|3T], [23 |. 

Let us explain now, following the reasonings of [31] and [23], that the 
masses of solitons are also generated by the spontaneous breakdown of the 
conformal symmetry. 



The energy momentum tensor of such theories is of the form (see |2TJ for 
more details) 

Lfxp = Q^u + Sfj, v . (95) 
where S^ v is the improvement term 

-V = -j Tr (Q s fa (B-'d, B) - 9ilv d p [B- X d p B 



I 

kN s 
I 



(d,d u - 9liu d 2 ) (96) 



Due to the fact we are dealing with a conformally invariant theory, L^ v 
satisfies 

d_ L ++ = , d+ = , L+_ = L_+ = . (97) 
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Even though it is not traceless, is symmetric and conserved, 

^0^ = 0, (98) 

The energy of classical solutions are given by the space integral of the 
(0,0) component of energy-momentum tensor L^ v . In the Lorentz frame 
where the classical soliton solution is static, the energy should be interpreted 
as the mass of the soliton. However, since the theory is conformally invariant, 
it has no mass scale, and the soliton mass should vanish. When the conformal 
symmetry is spontaneously broken by choosing a particular constant solution 
for the field rj, we obtain a massive theory. Construct the energy-momentum 
tensor of such a theory as follows. Clearly, the tensor M!/ , introduced in (|9~5|) 
and evaluated at any classical solution, satisfies (0). Therefore, the tensor 
defined by 

pvbroken _ a. I (qq\ 
^ txu — ^fiu |?7=constant j 

is symmetric and conserved, 

^Qbrokcn = q ^ ( 10Q ) 

since rj = constant is a solution of the equations of motion. Then, let the 
energy in the massive theory be proportional to the space integral of @Qo° kcn . 
Using ( p6] ) and (|95|), we obtain the soliton mass in the form 

Mr 00 k N 
^== = - dx 6^ oken + £ vac . = -^d x {v + Qx + x_) rj^Ol) 

V 1 — V J-oo I 

because the integral of Uqq' vanishes by the above arguments. Here v is 
the soliton velocity in the units of the speed of the light. Notice that we 
have subtracted the energy _E vac of the vacuum solution which is, in fact, 
divergent. Of course, the vacuum solution is not unique, and it is not clear 
which one provides the absolute minimum of the energy. We will use the 
following prescription for the soliton mass formula. For the soliton solutions 
lying, under the dressing transformations, on the orbit of the vacuum solution 
(|62"D, we take Q in ( |101| ) to be that one given in (|63|). However, for those 



soliton solutions lying on the orbit of the vacuum (|60|) , we take f2 in ( |101| ) 
to be equal to the parameter (3 introduced in floTil). Such a prescription 
guarantees the finiteness of the soliton masses. 
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The soliton masses are determined solely by the behaviour at x = ±00 
of the space derivative of the field v. That is quite a remarkable fact. In 
addition, as we now explain, it is very easy to obtain such a behaviour in the 
general case from the solitonic solutions fl87|). 

Consider an integral gradation of Q, with s[ = ^ s«, «o = — an d 

i 

Si labeling the gradation that defines the model (|55|)-([58]). Consider the 
integrable highest weight representation with highest weight state 

I K) = ® I A,) 8 *, (102) 

i=0 

where | A*) are the highest weight states of the fundamental representations 
of Q, and Aj are the corresponding fundamental weights of Q. 

Then it is possible to show |21[] that taking the expectation value of both 
sides of (^) in such state, one gets (with the gauge choice /il(x_) = }ir(x + ) = 
1) 

e -(u+nx +x .)N s 4 = ^ 1 e x + e+ e -x.e- pe x.e. e -x+e+ | Xgi y ( 1Q3 ) 

Now, choosing p to be the exponential of an eigenvector of £±, 

[S±,V}=u±V, (104) 

we obtain a soliton solution 

e -(, + nx + x.)N s ^ = J e e^V J Ag/) (1Q5) 

with T = u + x + — UJ-X- = 7 (x — vt). 

Suppose V is an operator in such a representation for which there is a 
positive integer N' Vl such that 

(As/ I V n I A s /) = for n > N' v . (106) 

Then the soliton mass is easily obtained from ( p.01| ), where for 7 > (7 < 0) 
only the upper (lower) limit x = 00 (x = — 00) contributes in the integral^ , 

2 k N' v . . , rr 2 2 k N' v , . . 



4 We point out that the soliton mass formula (107) could be equally obtained by defi ning 
the mass through the momentum formula, instead through the energy like in (|l0l| ), as 
= 2 = f da;0Q™ kcn . In this case, we do not have to subtract the vacuum momentum, 

since it vanishes. 
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Notice that we must have uj + uj- > in order to have the soliton velocity 
v = (c<j_ — u) + )/ (u- + u)+), not exceeding the light velocity (c = 1). 

The soliton mass formula (|107| ) has some remarkable properties. One 
of them concerns the relation particle-soliton in the theory, indicating some 
sort of duality similar to the electromagnetic duality of some four dimensional 



gauge theories possessing the Bogomolny (monopole) limit jn|. As we have 
seen, the soliton solutions are created by the eigenvectors V of S±. From (|104|) 
one has [£+ , [£_ , V]] = u + U-V. Expanding V over the eigenvectors of the 
grading operator Q s as V = Y,n V^ n \ one observes that [£+ , [£_ , V^]] = 
w + w_FW. Therefore, if some e Q n , n = 0, ±1, ±2, . . . ± (I - 1), does 
not vanish, it implies that must be one of the eigenvectors X in ( ]93|) . 
Then we associate a soliton with a fundamental particle. In addition, we 
have A = and, consequently, from (0) and ( |107|) , the masses of the 

corresponding soliton and fundamental particle are determined by the same 
eigenvalue. In fact, we have from (|93|) and ( |107| ), with v v = 1, that 

M„. = ^mr (108) 

Of course, in the expansion of V, we may have more than one non vanishing 
V^ n \ with n = 0, ±1, ±2, . . .±(1 — 1). Then we would associate a one-soliton 
solution to more than one fundamental particle. The counting of one-soliton 
solutions has to be better analysed in each particular case. We discuss this 



issue in section 10. 



9 The matter fields 

It is clear from (|65|)-(|69|), that the massive fields associated with non vanish- 
ing grade (namely F^), are chiral fields with non vanishing spins, in contrast 
with the Toda type fields. In fact, we show that the free equations for such 
fields take the form of the massive Dirac equation, as could be expected from 
general covariance arguments. 

Consider the subspace Q m for < m < I. Let Gm ' be the subspace of Q m , 
generated by the eigenvectors of [E_i , [Ei , •]] with non zero eigenvalues, i.e., 

= {T^ e Gm I ^ 0}, (109) 
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where A*" 1 -* is defined as 

[E^ , [E t , T^}} = [E t , [E-i , T™]] = TK (110) 
Decompose the subspace Q m , as a vector space, into the sum 

§ m =§yp+§gp, (in) 

where Qffi j s the complement of Qtfp in Q m . 

It is possible to show (see ^Tf for details) that the subspaces Gfi +m and 
Gffi are isomorphic. The mapping is given by the action of E-i on G^\ or 
equivalently by the action of E\ on Q_ l+m . Therefore, we can put in one-to- 
one correspondence the fields in G^-i+m and Qffl '. Then, one can show that 
each pair of such fields constitute a Dirac spinor under the two dimensional 
Lorentz group, and their equations of motion can indeed be written in the 
form of (obviously not free) Dirac equations. Consequently, we interpret the 
massive fields associated to generators of non zero grading as matter fields 



10 An example of a special class of models 

There is a class of models possessing a U(l) Noether current, which, under 
some circunstances, is proportional to a topological current. That occurs for 
those models where the grade I of the operator Ei, introduced in fl48f) , is 
equal to the integer N s defined in (0). In addition, it is necessary that the 
operators E± Ns satisfy the condition 

zE_ Ns = ^iz~ 1 E Nb G center of Go (H2) 

where, fi is some constant independent of z, and z is a complex variable used 
to realize the generators of the affine Kac-Moody algebra G, in terms of those 
of the finite simple Lie algebra G as 

H: = z n H a , E2 = z n E a , D = z^- (113) 

az 

That means that the "projections" of E±n s onto Go, are parallel and lie in the 
center of Go- When those condition are satisfied it is possible to gauge fix the 
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conformal symmetry, such that a special U(l) Noether charge is proportional 
to a topological charge. 

In this section we discuss a example where that happens. It correspond 
to the principal gradation of s/(2)^' with / = 2. Let us denote by H n , E±, 
D and C the Chevalley basis generators of the sl{2)^> . The commutation 
relations are 

[H m ,H n ] = 2mCS m+nfi , [E™, E n _] = H m+n + mC 5 m+nfi ,(lU) 
[H m , E±] = ±2E^ +n } [D,T m ]=mT m J T m = H m , E%; (115) 

all other commutation relations are trivial. The grading operator for the 
principal gradation (s = (1, 1)) is Q = \H° + 2D. Then the eigensubspaces 
are G = {H°,C,Q}, Q 2n+1 = {E™, E™ +1 }, with n e Z, and Q 2n = {H n }, 
with n G {Z - 0}. 

The mapping B is parametrised as 

B = e* H ° e* c e^ = e vC e^, (116) 

where H° = H° — ~ C is the Cartan generator in the special basis introduced 
in and so v = v — \tp. 
In the case I = 2, we choose 

E 2 = mH\ E^ 2 = mH-\ (117) 

where m is a constant. We then have 

[E 2 , [E_ 2) El]]=Am 2 El. (118) 

Therefore, each of the subspaces Q±\ has two generators with the same eigen- 
value 4m 2 . Following section |9| we write 

F+ = 2Vml (tPr El + $ R El) , Ff = 2yfml (iJ> L E+ 1 - ^ L E°) , (119) 

and introduce the Dirac fields 

*-(£)' *-($) (120) 

From (0) we obtain the masses of the particles, 

m v = ffij = = 0; = 4m; (121) 
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The equations of motion derived from (|55])-([58]), are 

1 



d\ = 


-4m^ 5 e v+2 ^ 5 ip, 


d 2 v = 




d 2 rj = 


o, 


vfd^ = 


m^p e v+2lp ^ip, 


vfd$ = 





(122) 
(123) 

(124) 
(125) 
(126) 



where the gamma matrices are defined as 



(i27) 

and 75 = 7o7i, and ip = V ;T 7o- Recall that d 2 = d 2 — d 2 , x± = t ± x. The 
corresponding Lagrangian has the form 



\c = ld^d^ + \d^d^+ l -d^d^- l -mle 2 ^ 

+ i^d^ -m^^e^ 2 *^^. (128) 



It is real (for 77 = real constant) if ip is the complex conjugate of ip, and if ip 
is pure imaginary. This will be true for the soliton solution as we shall see 
below. 

Notice that such model is invariant under the transformations 

x + <-> x- ; ip R <-> eV'L ; $r <-> -e^x ; <f ^ f, V 77 ; z/^z/ 

(129) 

where e = ±1. It should be interprated as the product CP of charge conju- 
gation times parity. Parity alone is clearly violated. 

The generator H° e Q commutes with E± 2 , and, therefore, the gauge 
symmetry ©-(^D of the model is U(l) L ® U(1) R , 

h L (x-) = e^ x - )H ° , h R (x+) = e Mx+)H °. (130) 

Since the genereators of U(1)l and U(1)r are the same, we have the global 
gauge transfomations ( [74]) generated by h D = h L = h]^ = e ldH °l 2 (9 = 
const.). The fields are transformed as 

ip — ► e l V ^ — > e~ 10 ^ (p — > (f , z>-^z>, 77 — ► 77; (131) 
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and the corresponding Noether current is 

J^ = <0 7 "<0, d„J^ = 0. (132) 

The fields if) and if) have charges 1 and —1, respectively; and tp, v and rj have 
charge zero. 

Let us next see how the general arguments given above concerning Noe- 
ther and topological charges apply here. The topological current and charges 
are 

f = ^-.e^d v V , Q topol . = [ dxf , (133) 
Zixi J 

Indeed, the Lagrangian ( |128| ) has infinitely many degenerate vacua due to 

the invariance under (p — > <p + iix. Making use of the field equations, one 

easily verifies that 

#757^ + o^^V 



9, 



y r (134) 

Combining this equation with the conservation of the vector current ifi^if), 
one deduces that there exist two charges defined by 

J = -i$RipR + ^9 + ip, J = ilpL^L + 

which satisfy d-J = 0, d + J = 0. We now make a "gauging fixing" of the 
conformal symmetry by choosing J = J = 0. We call it a "gauging fixing" , 
because any values of J and J can be transformed to zero by a conformal 
transformation. This gives, altogether, 

l-^ v d v ^=-i)^ifj, (135) 

ZTTt 71 

so that the topological and Noether currents are proportional. As discussed 
at the beginning of this section, that is a consequence of the fact that E± 2 
satisfies QTT2D. 

Let us turn to the Noether charge which here is simply the fermion num- 
ber. It should be defined such that it satisfies the Poisson bracket relation 

i {lf>, QNoetherjp.B. = ^ ( 13 6) 

Since the coupling constant k was taken as an overall factor, this is satisfied 
by 

^Noether = k J dxip>y ip (137) 
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so that | 

Qtopol. = ^QNoether (138) 

As argued in general, this means that k should only take discrete values as 
expected, since our actions are related with the one of WZNW. 

Let us now construct the soliton solutions. The operators E± 2 given in 
( 117 ), lie in the homogeneous Heisenberg subalgebra generated by H n , with 



the commutation relations (|114j) . Such a subalgebra has no generators of 
grade ±1 for the principal gradation. Therefore, the model under consider- 
ation has no vacuum solutions of type (Bp). Then, from (|76|), we get 



E± = E± 2 = mH ±1 . (139) 

We perform the dressing transformation starting from the vacuum solution 
(|60l), namely 

ip = r] = ip = ij=0, ; v = -imjx+x. = u . (140) 

Now, let | Ao ) and | Ai } be the highest weight states of two fundamental 
representations of the affine Kac-Moody algebra 31(2)^, respectively the 
scalar and spinor ones. Then, from (|84|) with rj = 0, we obtain the solutions 
on the orbit of the vacuum ( |140| ), 

e"* = (Al 1 ° 1 H = (Ao | G | A ), 
(Ao | G | Ao)' 

[m ( Ao I El G I Ao ) ~ fm ( Xi I E° + G I Ai ) 

Wr = \ -■ 71 — , -z , t ; » Wr 



(A | G | A ) V z (Ai | G | A x ) 

/m (Ai | G E°_ | Ai) ~ fm (Xo\GE^\Xo) 

where 

G^e^+^e-^^pe^^e-^^. (142) 

In order to get the soliton solutions, we choose the fixed mapping p to be 
an exponentiation of an eigenvector of E± (solitonic specialization); namely, 
p = e v , with [£± , V) — u± V. Therefore, 

G = exp (e r V^j with T = oo + x + — lu-X- = 7 (x — vt) . (143) 
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In this case the eigenvectors of £± are 



Indeed, 



i40)= E z ~ nE ±- ( 144 ) 

neZ 



[£ + ,V±(z)\ = ±2mzV±(z)=u$V±(z), (145) 
2tt? 

[£_,V±(z)] = ±—V±(z)=u£V±(z). (146) 
The solution, associated with V + (z), is 

!/ = u , v = $=0, ^ = ^fe r f _i J ; (147) 

while those, associated with VI (z), is given by 

m _ r ( 1 



u = uo, ip = ^ = 0, ^ = -J-e- L , , (148) 



where 



r = 2m(zx + x) = 7 (x — vt) . (149) 

z 

The masses of these solutions are obtained from ( |107| ). Here the relevant 
state | A s /) in (|106|) is 

I A B /> =| A )® | Ai). (150) 
Using level one vertex operators, one can verify that 

(A; | {V±(z)) n | Ai) = , for n > 1 and % = 0,1. (151) 

Therefore, Ny = in (|106|) , and from (|107|) one gets that the masses of the 
solutions ( |147| ) and ( |148[ ) vanish. Such solutions correspond to the objects 
which travel with velocities v — ± (1 — z 2 ) / (1 + z 2 ); and keeping z 2 > 0, 
one has | v |< 1. Therefore, these solutions cannot be interpreted as solitons 
(particles), since they would correspond to massless particles traveling with 
velocity smaller that light velocity. We should interpret them as vacuum 
configurations, since they have the same energy as vacuum ( 140 ). 

The true soliton solutions of the system are constructed as follows. Notice 
that V + (z) and V-(—z) have the same eigenvalues. Therefore, any linear 
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combination of them, leads to solutions traveling with a constant velocity 
without dispersion. So, we let 

V{a ± ,z) = Vi{a + V+{z) + a-V-{-z))] (152) 

[£ + ,V{a ± ,z)] = 2mzV{a ± ,z), [£_,V{a ± ,z)] = —V(a ± ,z), (153) 

and so uj + = 2mz and = — . The particular factor \fi is chosen such 
that the reality condition will be obeyed with a_ = a* + . Again, using level 
one vertex operators, one can verify thaif] 

(A s - | V(a±, z) n | A s /) = for n > 4. (154) 

Therefore, N v = 4 in (fLOg), and from ( W<\ ) with i/j 2 = 2, and ip being the 
highest root of s/(2), one gets that the mass of such solutions is 

M = 8km = 2km<p, (155) 

where k is the coupling constant appearing in the Lagrangian (|128|) . The 
solutions generated by (|1 52|) , have two parameters, namely a±. One param- 
eter is always present, because one can scale an eigenvector of £± without 
changing the width 7 and velocity v of the soliton, obtained from the eigen- 
vectors u±; see fl3T|). However, in this case, the second parameter comes 
from a symmetry As we have pointed out in (j39p, associated to the fixed 
element p = e v{a± ' z \ we have an orbit of equivalent solutions due to the 
global transformations ( |131| ), 

V(a±,z) -> Vi(a + e i9 V + {z) + a. e~ ie V.(-z)) . (156) 

The explicit form of the solutions generated by ( |152| ), is obtained using ( |141| ), 

tp = log -—; ), (157) 



tae 2r 



v = — \og(l + icre J — -m^x+X-, (158) 
T) = 0; (159) 



a Notice that the truncation occurs for powers greater than 4, and not 2, because | A s /) 
lies in the tensor product representation, see (150) 
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and 



ip = a. 



m e 



l—iae 2 



l—icre 2T 

-1 

l+iae 2r 



(160) 



where T is given in (|149|) , and a = a + a^z/A. Keeping m and z real, we 
have the mass M of the soliton, from ( |155| ), real and positive, and also the 
parameters 7 and v (|149| ) are real. The reality condition is obeyed if a_ = a+, 
as anticipated. At this point, it is useful to re-express the expressions just 
given in terms of the physical parameters of the soliton. Using equations 
(|127|) and (|121|) and one deduces that 



7 = I VI - v 2 , z = \J{1- v)/(l + v). 



Moreover, since a± are complex conjugate, we may write 



a± 



(161) 



(162) 



The dependence upon space-time appears only through y / aexp(r). We will 
write [] 

\fa^ = exp((7(x — x — vt)) 



where xo is the position of the soliton at time zero. Then we have 

ip = 2% arctan (exp (2m^ (x — x — vt) / y/l — v 2 j \ , 
which is the sine-Gordon soliton. The Dirac fields are given by 



(163) 
(164) 



J6 



im^ e ' 



; (x — XQ — vt) I \J\ — V 2 



( 



\ 



(l+l) 



1/4 



1+v 
1-v 



1/4 



(165) 



1—ie 



2m^(a; — XQ—vi)/ J yl—v^ ' J 



and ip is the complex conjugate of ip. Thus the only parameters are the 
soliton mass and velocity, together with the angle 6 which reflects the global 
invariance ( |131| ). Notice that the sign of the tolopogical charge can be re- 
versed by reversing the sign of z. Therefore, the solutions (|164|) - (|165 ) contain 
the sine-Gordon soliton and anti-soliton. 



3 by convention, we choose a to be positive 
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Finally, we come to the very important feature of the present model al- 
ready mentioned above in general, namely it is clear from the explicit expres- 
sions Eqs. |165| that if) vanishes exponentially when X — Xq — > ±00, so that the 
Dirac field is confined inside the soliton. That this must be true is of course 
a general consequence of Eq J135| which may be verified directly on the ex- 
plicit solution. This phenomenon has been much studied for electron phonon 
systems. Models of a similar type describe the electron self-localization in 
quasi-one-dimensional dielectrics (for recent reviews see p2 |, ||33|| ). At low 
temperature these systems go over to dielectric states characterized by charge 
density waves which can be constructed on the basis of the Peierls model. 
The continuous limits are described by Lagrangians similar to Eq. |128| . Dis- 



cussing this important issue is beyond the scope of the present article, so 
we will not dwell upon it here. Let us simply recall that the typical ex- 
ample of the polyacteline molecule was much discussed in connection with 
fermion number fractionization Clearly, on the other hand one may 

regard our soliton solution a sort of one dimensional bag model for QCD. 
In this connection let us note that, if we introduce the two- by-two matrix 
U = exp(r^ + 2^75), we may rewrite the Lagrangian Eq. |128| as 



C 



— (tr 
16 I 



mt 



-tr 



'U-%U] tr [U^d^u] } 

(166) 



which is similar to a two-dimensional version of the low energy effective action 
for QCD (see e.g. H). 



11 Conclusions 

From the common description based on zero curvature equations, a unified 
treatment has been given of the procedures to obtain solutions of nonlinear 
equations for a large class of integrable hierarchies including the general- 
ized KdV and non-abelian Toda coupled to matter. Interesting examples of 
them have been discussed in detail, taking advantage of the simplifications 
achieved, and the same procedure can be equally well applied to more com- 
plex equations. The role of the tau functions has been clarified, as specific 
matrix elements in special integrable highest weight representations (for any 
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level) . The results are a further step towards stablishing the conjecture that 
all multisoliton solutions lie in the orbit of the vacuum, generated by dress- 
ing transformations, which are now clearly related to the tau functions and 
solitonic specialization. These results should be very useful for the classifi- 
cation of integrable theories in two dimensions and for the generalization to 
the quantum case and to higher dimensions. 
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